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ABSTRACT 



Aims. We study the spectral properties of the unresolved cosmic X-ray background (CXRB) in the 1.5-7.0 keV energy band with the 
aim of providing an observational constraint on the statistical properties of those sources that are too faint to be individually probed. 
Methods. We made use of the Swift X-ray observation of the Chandra Deep Field South complemented by the Chandra data. 
Exploiting the lowest instrument background (Swift) together with the deepest observation ever performed (Chandra) we measured 
the unresolved emission at the deepest level and with the best accuracy available today. 

Results. We find that the unresolved CXRB emission can be modeled by a single power law with a very hard photon index r=0. 1+0.7 
and a flux of 5+^^x10"'- erg s ' cm deg-^ in the 2.0-10 keV energy band (Icr error). Thanks to the low instrument background of 
the Swift-XRT, we significantly improved the accuracy with respect to previous measurements. 

Conclusions. These results point towards a novel ingredient in AGN population synthesis models, namely a positive evolution of the 
Compton-thick AGN population from local Universe to high redshift. 



1. Introduction 
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The flux limit reached by the Chandra deep field South (CDF-S) 
4 Ms observation is such that the observe d source angular den - 
sity is four sources per square arcminute (iLehmer et al.ll2012l) . 
These sources have been found to be mostly unobscured and 
Compton-thin (Nh < lO^-' cm AGN at redshift < 2 with some 
contribution at soft energies (< 2 keV) from ga laxy clusters and 



starburst galaxi es (Brandt & Hasinger 2005; Tozzi et al] l2006t 
uo et al.ll200 8': Xue et al. 2011). Stifl, it is worth investigating 
whether a non-negligible contribution to the cosmic X-ray back- 
ground (CXRB) has not been identified yet. This may be due to 
an intrinsically diff'use component or to the cumulative contri- 
bution of individual sources below the current flux limit. In the 
latter case, high redshift and Compton-thick AGN (the ones for 
which the neutral hydrogen column density is higher than the 
inverse of the Thomson cross-section, Nh > 1.5x10^'*), together 
with star forming galaxies, are expected to be the main contrib- 
utors. 



The very fact that bright quasars, 
sive black holes, have been found 



powered 
z>6 



at 



supermas- 

(iFan et alJ l2006t 

IWillottetal] [2010: Mo rtlocket"al]|2011h . imphes the existence 
of a large numbers of less massive AGN at these early epochs 
dVolonteril llOlO). Thus, because no source at z > 6 has been de- 
tected in the CDF-S, the unresolved CXRB can be used to probe 
SMB H formation models at these early epochs (ISalvaterra et al.l 
l2012h . In the same way, the unresolved CXRB can be used to 
put some constraints on the contribution of the X-ray sources to 
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re-ionization jPiikstra et al.ll2004t ISalvaterra et"al]|2005L 120071: 
iMcOuinnllMa ! 

Although CDF-S observation reaches a flux limit such t hat star- 
burst galaxies match AGN in surface density ( Le hmer et al.l 
120 12h . most of the galaxies still remain undetected at low red- 
shift, contributing to the unresolved CX RB. Indeed, u sing the 
X-ray galaxy luminosity as proxy (Minec Tet alj|2012h . CXRB 
can be used to study the evolution of the star formation rate up 
to redshift z-6 (Dijksti-a et aL.,2012:,Cowie et al. 2012.) . 

Compton-thick (CT) are expected to be the main contributors 
to the unresolved X-ray emission in the CDF-S. In fact, while 
their n umber is expected to be similar to m oderately obscured 
AGN (iGilh et al.ll2007t iTreister et al.ll2009l) . only a few have 
been found among the 740 sources of the CDF-S (iTozzi et al.l 
2006; Luo et al. 2011) 

From an observational point of view two different approaches 
have been pursued to probe the unresolved emission in the CDF. 
The first approach consists in calculating the difference between 
the total CXRB and the i ntegrated flux of all detected sources. 
iMoretti et al] (l2003h and IWorslev et all (12006) . following this 
approach, found that the unresolved fraction was 10-20% , and 
higher at harder energies. Moreover. IWorslev et al.l (l2004l) rec- 
ognized the signature of a highly absorbed AGN in the spec- 
tral s hape of the unresolved emission. Unlike these, iLuo et all 
(1201 1) , while measuring similar values in the softer band, found 
that the unresolved emission is consistent with zero, in the 6-8 
keV band, posing a strict upper limit. The main sources of uncer- 
tainty here are the total CXRB measure CRevnivtsev et al. 20()5|; 
iFrontera et al.i 2007; M oretti et a n i2009h and the statistical con- 
tribution of the rare and bright sources. 
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The second approach is a direct me asure of the unresolved emis- 
sion :|Hickox^M^k^ito3<l2QQ7') found that the flux, with large 
uncertainties (~200%), was consistent with zero in the 2-5 keV 
band both in the CDF-S and CDF-N. The source of this huge 
uncertainty is the very low surface brightness of the unresolved 
emission, which represents only a small fraction of the Chandra 
instrument background. In this regime the signal/background ra- 
tio is so low that a few percentage points of systematic error in 
the background measure irreparably affect the measurement ac- 
curacy. Q 

In this paper we work out a different approach, directly measur- 
ing the unresolved emission on the Swift XRT observation of 
the CDF-S and subtracting the signal of the sources revealed by 
Chandra. In this way we exploit, at the same time, the low and 
predictable instrument background of the Swift X-ray telescope 
(XRT) and the unprecedented depth of the Chandra observation. 
The paper is organized as follows. In Section |2] we briefly de- 
scribe the Swift XRT telescope and the datasets we used, show- 
ing their distinctive qualities. In particular we dwell on the in- 
strument background comparing Swift XRT and Chandra in or- 
der to emphasize the peculiarity of our work. In Section [3] we 
deal with a wealth of technical details to give the full particulars 
of the spectral analysis used in this paper. For the sake of clarity 
the fit procedure and the results are presented in Sect. |4] Finally, 
in Section |5] we discuss our results, first, comparing them with 
previous measurements and then, with the expectations from one 
AGN population synthesis model. 

Throughout this paper, errors are quoted at the 68% confidence 
level for the two parameters of interest (A;t'^=2.3), unless oth- 
erwise specified. We adopt the following cosmology parameter 
values: Q.^ = 0.3, Qa = 0.7, ho = 0.7. 

2. Instrument, data, and reduction procedures 

2.1. Swift-XRT 

The X-ray telescop e (XRT) on board the Swift satellite 
jGehrels et al.ll2004 . uses a Wol ter I mirror set, or iginally de- 
signed for the JET-X telescope jCitterio et al.l 11994 ). to focus 
X-rays (0.2-1 keV) onto a XMM-Newton/EPIC MOS CCD de- 
tector dBurro ws et al. 2005). The effective area of the telescope 
(~ 150 cm^ at 1.5 keV) is about five times smaller than Chandra. 
The PSF, similar to XMM-Newton, is c haracterized by a ha lf en- 
ergy width (HEW) of ~ 18"at 1.5 keV (iMoretti et al.ll2007l) . The 
nominal field of view has a radius of ~ 1 1 arc minutes, with a 
pixel scale of 2.37" 

2.2. XRT observation of the CDFS 

Swift observed the CDFS for a total nominal exposure of 563700 
seconds in the period July 2007 to December 2007. We reduced 
data using the standard software (HEADAS software, v6.10, 
CALDB version Sep 10) and basically following the procedures 
described in the instrument user guide with some minor mod- 
ifications for our purposes. In our analysis we excluded the ex- 
ternal (detx < 90 and detx >510 ) CCD columns that are af- 
fected by the presence of out-of-time-events from corner cal- 

' This is exactly the same situation as was faced in the study of the ex- 
ternal regions of galaxy clusters where the intergalactic medium (ICM) 
surface brightness is at the same level. In this kind of study telescopes 
such as Suzaku and Swift are preferred to Chandra and XMM, which 
have much high er instrument background and a less favora ble ratio sig- 
nal/background dBautz et al.ll2009l : lEttori & Molendill2010l) . 

^ http://heasarc.nasa.gov/docs/swift/analysis/documentation 
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Fig. 1. The energy channel (PI) distribution of the XRT unre- 
solved emission (black) compared with the instrument back- 
ground (red). In this plot PI channels have been transformed in 
energy using a single value and not the RMF matrix. 



ibration sources. This corresponds to a nominal field of view 
(FOV) of 16.5 X 18.9 arcmin (0.087 deg^). To detect sources 
we ran the task wavedetect (Freeman et al. 2002) with default 
input parameter values on three different energy bands (0.5 - 
2.0; 2.0-7.0; 0.5-7.0 keV), on a large number of spatial scales 
(2,4,8,12,16,24,32,64,128 pixels, 1 pixel equals 2.36"). Merging 
the three catalogs yields a total of 109 sources down to a limit 
of ~ 3 xlO"'^ erg s"' cm"- and ~ 1x10"'^ erg s"' cm"^ in the 
0.5-2.0 and 2.-10. keV bands respectively. 

2.3. XRT instrument background 

For our purposes, the most relevant characteristic of the XRT is 
the low level and high reproducibility of the instrument back- 
ground (NXB), which we define as the signal registered on the 
CCD when focused radiation is excluded. This is mostly con- 
tributed by the in duced particle em ission and the pure instrument 
electronic noise (iHall et alJl2008l) . 

To estimate the NXB spectrum in the CDF-S observation, we 
used the signal registered on the small regions in the four CCD 
corners that are not exposed to the sky (NFS regions) and not 
contaminated by calibration sources. As these data are only 
available for observations after June 2008, we used observations 
performed in 2008-2009 for a to tal of ~ 7 Ms e c. For a deta iled 
study of the XRT NXB also see iMoretti et al.1 (l2009i 12011 . To 
check at what level of accuracy we can reproduce the XRT in- 
strument background we compared this spectrum with the XRT 
CDF-S observation in the 7-10 keV energy band. There the ef- 
fective area is a factor five to ten times lower than the one at en- 
ergies <2keV and, at the same time, the instrument background 
is inflated by the presence of the Ni (Ka and at 7.478 and 
8.265 KeV) and Au {ha at 9.713 keV) fluorescence lines. Once 
the detectable sources have been removed, the cosmic signal, is 
expected to be only < 1 % of the background in this particular 
energy band. 

As shown in FiglUwe find that the two spectra coincide very 
well beyond 7 keV. Indeed, the ratio of the signals between 7-10 
keV is XRT-CDF-S/NES = 1.033, while if we restrict to 7 and 8 
keV, where the nickel K lines are completely overwhelming, the 
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Fig. 2. Comparison between Chandra and Swift-XRT instrument 
background spectra. Upper panel: NXB of the two instruments 
are plotted per unit area, to show the influence of different envi- 
ronments (orbits plus satellite). Middle panel: NXB are plotted 
per solid angle, to show the influence of different focal lengths. 
Lower panel: the ratio between Chandra and XRT NXB calcu- 
lated per solid angle and normalized by the respective effective 
areas. This shows the improvement in the signal/background ra- 
tio for an extended source using XRT instead of Chandra. 



agreement is even better XRT-CDF-S/NES = 0.993. As such, 
in the systematic error calculation we conservatively assumed 
that the instrument background is known at 5% accuracy (Icr). 
We note that , in some previous works t hat used Chandra and 
XMM data ( De Luca & Molendil l2004t | Hickox & MarkevitchI 
I2OO6) , the ratio of the signals at the very end of the energy 
band has been used to rescale the measurement of the instrument 
background. Here we find very good agreement between the 
two datasets, without any ad-hoc renormalization. This plainly 
shows the high level of reproducibility of the XRT instrument 
background. 

2.4. Comparison with Chandra 

To compare the Swift-XRT and Chandra instrument back- 
grounds, we considered the signal registered on the XRT NES 
regions and the signal registered on the ACIS-I in stowed po- 
sitioi£l In the upper panel of Fig. |2] we plot the two instru- 
ment background spectra normalized per unit of area, which 
show the effect of the different environments (orbit plus satel- 
lite). The very fact of being in high orbit (assuming negligible 
effects from the satellite structure differences) produces a par- 
ticle background in Chandra that is a factor ~5 higher than in 
XRT, which is in a low orbit. In the spectroscopic study of ex- 
tended sources, what is relevant is the background per unit of 
solid angle (middle panel of Fig. |2]i, which can be calculated 



^ File acis-D-012367-stowed-psu-evt-041104.fits downloaded from 
|http://cxc.harvard.edu/contrib/maxim/acisbg/ 



from the previous one, by accounting for the pixel scale. Since 
the ratio of the pixel linear scale between Chandra and Swift is ~ 
4.5, we find that the Chandra background per unit of solid angle 
is ~ 50 higher than Swift XRT. Finally, if we take the effective 
area into account (which is approximately a factor 5 higher in 
Chandra), we find that the ratio signal/background for extended 
sources (for which the angular resolution is not important) is a 
factor ~10 better in Swift-XRT than in Chandra (lower panel of 
Fig.H. 



3. Spectral analysis 

We restricted our analysis to the sky circular region centered 
on ra=53.1092 and dec=-27.82086 with a radius of 0.1053 de- 
grees (0.03184 deg^), which optimizes the intersection of the 
two instrument exposure maps. Hereafter we refer to this re- 
gion as the region of interest (ROI) region. The average effec- 
tive exposures within this region is 508 ksec . Moreover, we 
restricted our analysis to the 1.5-7.0 keV energy band. We ex- 
cluded the soft part (<1.5 keV) because in this range our data 
did not allow us to disentangle the genuine extragalactic CXRB 
components from the local ones. In fact, in this energy band, 
the diffuse X-ray emission is expected to be contributed both 
by the thermal e missi on from the Galaxy an d the local Bubble 
dSnowden et"al]|1998l:lKuntz & Snowdenli2()0C i) and by the ther- 
mal emission from faint galaxy groups and WHIM filaments 
(Cappelluti et al. 2012; Shull et al. 2011). Above 7 keV, as said 
in the previous section, the instrument background and, in par- 
ticular, the particle-induced one overwhelms the cosmic signal. 
To study the spectrum of the unresolved CXRB on the Swift- 
XRT data, first, we excised the signal of all the detected sources, 
which are 37 within the ROI, with fluxes in the range [le-15,le- 
14] erg s ' cm"^ in the 2.0-10.0 keV band. To get rid of most 
of the source signal without losing too much area, we excluded 
a circular region with a radius that depends on source counts. 
We used a radius equal to twice the radius at which the source 
surface brightness profile matches the background: this means a 
typical radius of 25 arcsec with values ranging from 20 to 60. 
Hereafter we refer to this radius as Xext ■ 

Besides the instrument background, the unresolved signal on the 
XRT detector in the 1.5-7.0 energy band is expected to be con- 
tributed by the following elements: 

- the sources detected by Chandra and not by XRT, 

- the optically/IR detected (X-ray undetected) sources, 

- PSF residuals from detected sources, 

- the stray-light contamination, 

- the unresolved CXRB that is the goal of the present work. 

To perform spectral analysis we subtracted the instrument back- 
ground and modeled the remaining five components. In the fol- 
lowing we describe item by item the way we quantified each 
single component. 

3.1. Chandra sources 

We used the public catalog of IXue et al.l (1201 ih . which is the 
result of the cross-correlation of three different energy bands 
and which consists of 740 sources down to the unprecedented 
flux limits of 10"'^ erg s"' cm"^ and 10""' erg s"' cm"^ in 
the 0.5-2.0 and 2.0-10.0 keV band. In the selected region the 
Chandra catalog lists 326 sources that have not been detected 
by XRT. We summed the pha files of the sources and of the 
(locally extracted) backgrounds, weighting the response (RMF) 
and ancillary (ARF) files by the source counts. The full de- 
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Table 1. Spectral analysis results of 32 stacked XRT sources and 
corresponding Chandra sources. 





slope fluxl-2 keV[erg s ' cm 


Chandra 
XRT 


1.09±(0.01,0.01) 2.59+(0.03,0.03)E-14 
1.16±(0.07,0.08) 2.71±(0.19,0.18)E-14 



-27.65 




-27.95 r , , I , , , , I , , , , , , I , -1^ , , I 

53.00 53.05 53.10 53.15 53.20 53.25 
R.A.(J2000) 



50 years after the discovery. 




Fig. 4. The PSF model of the 109 sources detected in the XRT 
observation at 2.5 keV energy. The large white circle shows the 
boundaries of our ROI, while small black circles show regions 
we excised (rext)- The PSF residuals is the small, but not negligi- 
ble signal outside small circles and within the ROI. 



Fig. 3. Black dots are the 18296 sources from iGrazian et af] 
(12006) catalog. The green circle indicates our ROI. Red dots are 
the 5683 optical/IR sources we used in the spectral analysis. 

ta ils of the reduction of th e X-ray spectral data are described 
inlVattakunn el et al.l (l20I2h . 

To quantify the impact of cross-calibration between Chandra and 
Swift XRT and source variability (Chandra and Swift XRT ob- 
servation were performed in different periods), we compared the 
summed spectra of the sources detected by both instruments, 
excluding the five brightest. These are 32 sources with fluxes 
ranging between 10"'^ and 10"''* erg s"' cm"^ in the [0.5-2.0] 
band. The spectral analysis of the two stacked spectra yielded 
results that are consistent within Icr at the level of 7% (TablT]). 
This variance can be easily encompassed in the systematic error 
budget of the analysis. This should be considered as an upper 
limit as the total flux variation associated to the intrinsic vari- 
ability of the 294 remaining fainter source sample is expected to 
be significantly lower than the 32 brighter ones. At the end we 
conservatively adopted a 7% (Icr) value in the systematic error 
uncertainties in the cross-calibration between the two telescopes. 



3.2. Optical/IR galaxies 

Several studies put in evidence that optical/IR galaxies, which 
are not detected in X-ray, show significant emission once 
their contribution is averaged (o r stacked ) over a large sam- 
ple (Worslev et al. 2006; Hickox & Markevitchll2007l) . We used 
the iGrazian et al.. (2006) multi-wavelength (from 0.3 to 8.0 
fim) catalog, which covers the large and deep area in the 
GOODS Southern Field by combining public data from Hubble 
VLT, Spitzer, and 2.2ESO. We used the 5683 galaxies detected 
within the four arcminute radius circular region centered on 
R. A. =53. 126 and Dec=-27.820, off'set by 1 arcminute to the west 



from the center of the ROI (Fig. O. This ensures a coverage of 
40% of our ROI. We found that these galaxies make a low, but 
significant contribution to the X ray emission, with a ratio sig- 
nal/background of ~3%. To test the systematic uncertainty of 
this operation we performed the same analysis by varying the 
region (from 3 to 6 arc minutes of radius) and the extraction ra- 
dius for the single sources (from 2 to 4 pixels). We found that 
scatter of the result is 10% at Icr level, which is much less than 
the statistical error. 

3.3. PSF residuals 

As said, from XRT data we excluded most of the signal of the de- 
tected sources by excising circular regions, with the size depend- 
ing on the source flux (rext). This leaves a small fraction (~ 5%) 
of source fluence spread on the rest of the field o f view, which w e 
evaluated using the PSF analytical model ( Moretti et al.ll2005h . 
We also considered the small contribution from PSF wings from 
sources detected outside the ROI (76 objects). This fraction de- 
pends on the energy since harder photons are more scattered. We 
calculated the residuals on a grid of 50 different energies in the 
2.0-7.0 keV range. We found that, at 2.5 (6.0) keV, 3.7% (4.0%) 
of the 109 source flux is diffiise all over the ROI (Fig.Hl). To ac- 
count for this energy dependence we built an ad hoc ARF file, 
linearly interpolating the grid test energies, such that fitting the 
stacked spectrum of the 109 sources yields the appropriate cor- 
rections. In the systematic error calculation we assumed that the 
PSF model extrapolated up t o ~6 arcmin radii is accurate at the 
level of 10% accuracy (Icr), dMoretti et al.ll2005h . 

3.4. Stray-light contamination 

One of the main components in the unresolved signal in the 
Swift-XRT data is the stray-light contamination. This is pro- 
duced by photons coming from sources that are outside the tele- 
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dist. from cenfer [i 



Fig. 5. Upper-left panel: Fraction of the flux of a source outside 
of the field of view falling within the annulus between 5 and 6 arc 
minutes. Black, red, and blue points represent the fraction mea- 
sured in the Crab calibration observations in the total (1.5-7.0), 
soft (1.0-2.0), and hard bands respectively. Lines are the linear 
interpolation between the measures. Upper-right panel: Black 
line is the calibration measure eiTor histogram. Green dashed 
line is the histogram of the relative difference, between soft and 
hard band. Lower-left panel: The stray-light contamination on 
the XRT detector from a uniform source in different concentric 
annuli on the detector, accounting for the CDFS exposure map. 
Right-lower panel: Black line is the error histogram. Green 
dashed line is again the histogram of the soft/hard relative differ- 
ence. 



scope FOV at distances between 25 and 100 arcminutes from the 
optical axis of the telescope. A fraction of the photons produced 
by these sources reach the detector after only one reflection on 
the mirror or even directly, passing through the mirror shells 
without any interaction. Some X-ray telescopes mount baffles on 
top of the miiTors to prevent this contamination. This is not the 
case for XRT, for which the stray light is a significant fraction 
of the diffuse radiation registered on the CCD. This conta mina- 
tion has been accurately calibrated in lMoretti et al.l (l2009l) using 
a set of observations of bright sources (Crab, Sco-Xl) observed 
outside the field of view. 

To calculate the stray-light contamination in the XRT CDF-S 
observation, we first split the detector into 12 concentric annuli 
(external - internal radius - 1') and calculated which fraction 
of the flux of the calibration source, positioned at a given angle 
outside the field of view, falls within the ith annulus (hereafter 
fraz,(0)). For example, we found that, when the Crab is posi- 
tioned at 45' from the telescope aim, (0.25+0.01)% of its flux 
falls within the 5-6'annulus (upper-left panel of Fig.|5]). The typ- 
ical 1 cr relative eiTor is ~5% (upper-right panel of Fig. |5]l. As 
shown in the same figure we do not detect any dependence of the 
stray-light fraction on energy: the relative differences between 
the soft (1.0-2.0 keV) and hard bands (2.0-7.0) keV, defined as 
abs(fraz.5o//-fi"az/,fl,-rf)/fraz5o/,, are at the level of the statistical er- 
ror (-5%). 

Then, for each annulus we calculated the stray-light contamina- 
tion produced by a uniform source extended all around the field 
of view. To do this we linearly interpolated the calibration obser- 
vations in the 15'< 9 <120'interval and integrated fraz,(0) In 9 
d9 (lower-left panel of Fig.|5j- In this way we get the contami- 
nation measure in terms of CXRB fraction. The histogram of the 
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Fig. 6. The five data sets we used in this work. The XRT un- 
resolved spectrum (black) has been fit by a sum of 5 different 
power law; (i) Chandra sources, not detected by XRT (blue); (ii) 
X-ray emission of optical/IR sources (violet) ; (iii) XRT sources 
to model PSF residuals (red); (iv) the CXRB spectrum to model 
stray-light contamination (green). The XRT unresolved spec- 
trum has been offset up by a factor 8 for clarity. 



relative eiTors of the 12 bins are plotted in the lower right hand 
panel: typical value is 5%. 

After summing different contributions at different off-axis angles 
and accounting for the effective exposure map of our observa- 
tion, we found that in our ROI the stray-light contamination is 
expected to be equivalent to the 0.690+0.020% of the CXRB or 
(4.9+0.2) xl0"'2 erg s"' cm^^ deg^^ in the 2.0-lO.OkeV. This is 
done in the assumption that the CXRB all around the CDFS is 
uniform and accurately reproducible. The first assumption is jus- 
tified by the fact that the master ROSAT catalogs Q do not report 
bright sources within a two degree radius. This means that we 
can safely assume that the contamination on our field of view is 
produced by a typical collection of CXRB sources that are sym- 
metrically distributed. It also follows that the second assumption 
is fully justified. In fact, although we do not have any direct X- 
ray observation of the nearby regions (ROSAT observation are 
limited to energies < 2keV), the cosmic variance is < 3% be- 
cause the interested sky area is as large as -3 deg^ (iMoretti et al.l 
20_09l). These uncertainties are accounted for in the systematics 
calculation, see Sect l4.2l . 



4. Fit procedure and results 

We used a total of five different datasets (Fig.|6]l, plus one for the 
instrument background: (i) the Chandra stacked spectrum of the 
326 sources undetected by Swift-XRT; (ii) the Chandra stacked 
spectrum of the 5683 optical/IR sources; (iii) the XRT stacked 
spectrum of the 109 detected sources, used with the ad-hoc mod- 
ified response file to model the PSF residuals; (iv) the CXRB 
spectrum (stacking of 130 XRT observations) used to quantify 
the stray-light contamination; (v) the fifth dataset is the XRT un- 
resolved spectrum that we assume is given by the sum of the 
four other (opportunely renormalized), plus the truly unresolved 
X-ray emission. To correctly propagate the statistical uncertain- 
ties of the single components, we fitted the five datasets together. 



http : //www . xray . mpe . mpg . de/ ~ j er/rosat/cats/catcat . html 
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Fig. 7. Results of the fit procedure. We modeled the XRT un- 
resolved spectrum as the sum of five components, each fit 
by a power law. Colors are the same as in the previous fig- 
ure: (i) Chandra sources not detected by XRT (blue); x (iii) 
Bright sources PSF residuals (red); (iv) stray-light contamina- 
tion (green); (iv) X-ray emission of optical/IR sources (violet) . 
The fifth component, the truly unresolved emission, which is our 
aim, is plotted with its uncertainties as the gray area. 
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Fig. 8. The photon index of the unresolved CXRB as function of 
the flux of the excluded sources. 

Table 2. Total flux of the 326 Chandra sources not detected in the 
XRT image: comp arison of our bes t fit to the stacked spectrum 
with the sum of the lXue et al.l (1201 ll) catalog photometric values. 





fluXo.5_2.OkcV fluX2.o_8.OkeV 

[erg s"' cm"^] [erg s"' cm"^] 


Present work (spec.) 
Xue et al. (2011") (photo.) 


2.69t;j^xl0-i^ 1.12tJJ^xlO-'^ 
2.70+«»5xl0-'^ 1.05!»!!:5xlO-'3 



We modeled (i)-(iv) datasets with single absorbed power laws. 
For the dataset (v) we used the sum of five different power laws 
(wabs*(pow-Hpow+pow+pow4-pow) in XSPEC syntax): the pa- 
rameters of the first four power-laws were linked to the best val- 
ues of the (i)-(iv) dataset fit, while the fifth photon index and 
normalization were left free to vary. The absorbing column was 
frozen to the Galacti c value derived from the HI Galaxy map 
(lKalberlaetalJl2005h 7 .02x10 cm . In the fit procedure we 
used C-statistics grouping the spectrum to have a minimum of 
one count for each bin. We checked that enhancing the minimum 
event number for each bin (we tested 3 and 5) does not signifi- 
cantly affect the fit results. We also checked that grouping data 
with a minimum of 20 counts for each bin and using chi-square 
statistics does not significantly change the results. 

4.1. Results 

Considering the five data-sets together, the adopted model pro- 
vides a good description of the data: chi squared = 2260.25 using 
1954 PHA bins; reduced chi-squared =1.16 for 1944 degrees of 
freedom. The best values are reported in Tab. [3] and the best- 
fit models are plotted in Fig. Q Instrument background and the 
stray-light contamination represent ~ 52% and 15% of the sig- 
nal respectively. The remaining 33% is essentially contributed 
by the emission of the Chandra sources and the real unresolved 
emission at the same level (15% each), while the optical sources 
represent 2% of the total. 

We find that the unresolved emission can be modeled by a hard 
power law with photon index r=0. 1+0.7 and a normalization of 
7.6 (h-10.0,-5.2) xlO""** photons s"' cm"^ deg"^ at 1 keV, cor- 
responding to a flux density of 5 xlO"'^ erg s"' cm"^ deg"^ in 
the 2.0-10. keV band. We note that excluding fainter and fainter 
sources makes the unresolved spectrum harder. As shown in 
Fig. [8j considering the whole signal from the CDF-S, includ- 
ing all the sources, yields a photon index is 1.20+0.05, which 



is close to flie total CXRB (r=1.4). Excluding all the XRT de- 
tected source makes the slope 0.84+ 0.1, while excluding all the 
Chandra 4 Ms catalog sources brings the photon index down 
to 0.7 + 0.3. Finally removing the contribution of optical/IR 
sources makes the final r=0. 1+0.7. 

In the remainder of this section we detail the results we obtain 
for the four other components of the XRT unresolved signal. 

(i) We find that the integrated emission of the 326 Chandra 
sources is well constrained both in slope (0.94^q qj) and flux 
(4.94+0 ]2 xl0"'2 ei-g ^-i ^m-^ in the 2-10 keV band). As ex- 
pected, the total integrated emission of these faint sources ( 
flux < 10"'^ and 10"'^ erg s"' cm ^ in the 0.5-2.0 and 2-10 
keV bands respectively) is harder than the spectrum of brighter 
source s. As a check we compared these fluxes with the photom- 
etry of IXue et al ] (120111) catalog, which is reported in the 0.5-2.0 
keV and 2.0-8.0 keV bands. We find very good agreement be- 
tween the two measures both in the soft and hard bands (Tab.|2]i. 

(ii) Considering the sources detected in XRT to calculate the PSF 
residuals, we get a photon index of 1.61;^Q Qg. 

(iii) The best-fit slope of the 130 XRT stacked observations (to 
quantify the stray-light contamination) is 1.46 ^0 ^4. 

(iv) The X-ray emission of the optical/IR sources can be mod- 
eled well by a relatively soft power law (photon index 1 .9) and 
a flux of 7.2+2-^ xlO"'^ erg s"' cm^^ deg^^ in the 2-10 keV 
band. These values are in very good agreement with Cow ie et aP 

who use the same optical/IR catalog and report a flux of 
~ 7.0 xlO"'3 erg s"' cm'^ deg"^ in the 2.0-10.0 keV band, (cal- 
culated from their Tab. 1) and a photon index r= 1.7-2.0 based 
on X-ray colors (Fig.|9]l. Our measure is slightly lower butcon- 
sistent with lHickox & MarkevitchI (120071) and lLuo et all (l201lh 
(up to 6 keV) who use different catalogs. In the hard part of 
the spectrum, we cannot confirm the results of iLuo et"al] (1201 Ih 
who find a strong signal equal to 4.8±2.0 erg s"' cm"^ deg"^ 
(28%±10.0 of the total CXRB) in the 6-8 keV band. While their 
measure is in good agreement with our spectroscopy in the 0.5- 
6.0 band, consistent with a photon index r=2, their photometry 
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Fig. 9. Comparison among different measures of the contribution 
of optical/IR (X-ray undetected) to the X-ray background. The 
gray area represents the best fit to our data. Green,_blue, and 
red po i nts represent th e ph otometry of Hickox & MarkevitchI 
(l2007h . lLuo et alj(l201 ll) . and lCowie et al.l(l2012l) . Black squares 
and arrows represent our photometric values and 2cr upper lim- 
its including systematic error of 2% on the background measure. 
Gray arrows show the upper limits accounting only for statistical 




Energy [keV] 

Fig. 10. The spectrum of the unresolved CXRB: comparison 
with previous measures. The gray area is the result of the p resent 
work. Green p oints are fro m Hickox & MarkevitchI (l2007h : blue 
points are from lLuo et al.l (1201 Ih . wi th the blue arrow indi cating 
the Icr upper limit; red points from ' Wors lev et afl (12006 ). The 
dashe d area represents the total CXRB emission as measured by 
Swift dMoretti et al.ll2009l) . 



is an order of magnitude higher than our spectrum in the 6-8 
keV interval. They interpret this strong feature as the signature 
of a substantial population of highly absorbed AGN (IXue et al.l 
120121) . The discrepancy with our data cannot be explained by 
the d ifferent optical/IR adopted catalogs: indeed with their cat- 
alog (iGiavalisco et al.l 120041) we find results that are fully con- 
sistent with our previous finding both in spectral slope and in 
flux. In order to go into tho discrepancy thoroughly, we also per- 
formed the photometry in the [0.5-2.0], [2.0-4.0], [4.0-6.0], and 
[6.0-8.0] keV bandsO While we find agreement with all pre- 
vious works up to 6 keV, we do not find any significant flux 
beyond 4 keV (Fig|9]l- Indeed, allowing only for statistical un- 
certainties our photometry in the 6-8 keV band yield s an upper 
limit that is inconsistent with the flux measured by iLuo et all 
([2011). However, this discrepancy can be explained by a slight 
difference in background estimate. Indeed allowing for only 
a 2% systematic error in our background value, we find that 
our 2cr upper limi t beco mes consistent with the photometric 
point by iLuo et all (1201 ih . This is because, once all the bright 
sources have been removed, in the 6-8 keV, where a strong in- 
strument feature is also present, the ratio between the instrument 
background and the cosmic signal is < 1% (see for example 
^http://cxc. harvard.edu/contrib/maxim/stowed/). 



^ We calculated the signal registered within the 1" radius circular 
regions centered on th e positions of the 8637 sources listed by the 
IGiavalisco et alj l l2004l) catalog within the 4. 2' radius circular region 
shown in Fig.jS] in the [0.5-2.0], [2.0-4.0], [4.0-6.0], [6.0-8.0] keV bands. 
To estimate the background we use the mean flux recorded in the same 
region after excluding of both the X-ray and optical/IR sources. In the 
6-8 keV we count 25,385 total events with 25,331 expected background 
events over a total of 103973 pixels. To convert these counts in flux, we 
use a mean exposure time of 3296249.2 seconds and a conversion factor 
of 1 count s"' = 7. 1x10" "erg s"' cm"-. Using a larger aperture for the 
source extraction regions we find consistent results. 



4.2. Systematic errors 

Because we work in a very low signal/background regime (15%) 
in the uncertainties calculation we considered the impact of the 
systematics of the single components of the model. To do this we 
repeat the fit procedure 1,000 times, each time varying the nor- 
malization of the different components according to Gaussian 
distributions with the variances reported above: 5% for the in- 
strument background, 7% for the Chandra source spectrum, 5% 
for the stray light, and 10 % for PSF residuals. We found that 
the impact of these systematics on the final slope uncertainty 
is modest, 0.08 (Icr), which is only 10% of the statistical error. 
The impact on the flux uncertainty is greater at 1 .77 erg s"' cm"^ 
deg"^, which is 75% of the statistical error. 



5. Discussion 

5. 1 . Comparison with previous works 

The unr esolved emission has been di rectly measured at faintest 
level bv lHickox & MarkevitchI (l2007l) who found that in both the 
CDF-S and CDF-N it was consistent with zero, although with 
large uncertainties: in the CDF-S they measured 1.4+3.9 xlO"'^ 
erg s"' cm"^ deg"^ in the band 2-8 keV. This is nearly a factor 
three smaller than our measurement, but still consistent within 
Icr (Fig. [Tol l. One could find surprising it that the level of un- 
resolved emission we measure in the 4Ms is three times higher 
than the one found by a previous work based on 2Ms. The rea- 
son for this is that the sources in the 4 Ms observation that are 
not detected in 2 Ms supply only a small fraction of the total 
CXRB flux. Indeed, whi le there are many sources of the 4 Ms 
catalog dXue et alj|201 Ih that are not included in the 2Ms cat- 
alog bv Luo et al.f (1200 8*) (169 out of 430 in our ROI, that is 
40%), their total flux is only ~10"'"* erg s"' cm"^, which repre- 
sents only 1.5% of the total CXRB or ~ 7% of the unresolved. 
However we stress that our measure represents an improvement 
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slope 


flux [erg s ^ cm ~ deg " ] flux [erg s ^ cm deg ^ ] 
(1.0-2.0) keV (2.0-10.0) keV 


fract.[%] 
(1.5-7.0)keV 


Chandra sou. 
PSF res. 
S.L. contam. 
OPT/IR sou. 


0.94 (-0.03,4-0.03) 
1.61 (-0.05,+0.05) 
1.46 (-0.04,+0.05) 
1.90 (-0.37,+0.45) 


5.70e-13 (-1.39e-14,+6.63e-15) 4.94e-12 (-1.20e-13,-l-1.23e-13 ) 
1.84e-13 (-7.53e-15,+5.74e-15) 6.88e-13 (-2.81e-14,-l-2.90e-14 ) 
1.08e-12 (-3.95e-14,+3.00e-14) 4.79e-12 (-1.75e-13,+ 1.98e-13 ) 
2.76e-13 (-7.18e-14,+8.94e-14) 7.25e-13 (-1.88e-13,+2.58e-13 ) 


14.7 
2.0 
14.3 
2.2 


Unresolved 


0.09 (-0.78,+0.70) 


1.77e-13 (-1.16e-13,+ 1.51e-13) 5.00e-12 (-2.62e-12,+3.22e-12 ) 


14.9 



Table 3. Best fit values of the different components contributing the unresolved signal in the XRT observation. Errors calculated 
at 68% confidence for two parameters (l^x^-2.y). In the last line we report the power-law best-fit values to the extragalactic unre- 
solved background. The reported errors are the quadratic sum of statistic and systematic error The last column reports the relative 
contribution of the single components, taking into account that the instrument background is 53% of the total. 



in the sense that the Icr upper limit (measured value -H Icr error) 
is lower than the previous ones. 

In the soft band we find that the spectral fit between 1.5 and 2.0 
keVis only marginally consistent with the lHickox & MarkevitchI 
('2007') photometry in the band 0.5 -2.0. However, we note that 
a fair comparison between the two measures is not possible be- 
cause the energy bands do not completely coincide. Moreover, 
in the 0.5-1.5 keV band the spectrum of the unresolved back- 
ground is expected to have a thermal component owing to un- 
detected galaxy groups and WHIM emission. This means that in 
the soft band the extrapolation of our best fit probably underes- 
timates the CDFS unresolved background. As already said we 
limit our measure to energies higher than 1 .5 keV because our 
data do not allow us to accurately disentangle these thermal ex- 
tragalactic components from the local thermal emission of the 
Galaxy and Local Bubble. 

In Fig. [TOl we al so plot the results of IWorslev et al.l (l2006h and 
iLuoet 311(1201 Ih . which are in good agreement with our results 
in the 2-6 keV band. For the reasons already discussed our spec- 
trum is not directly comparable to their photometry in the 0.5 
-2. band. 

iLuo et al.l (1201 Ih find that in the 6-8 keV the unresolved 
emission is consistent with zero with an upper limit at 10% of 
the total CXRB at 7 keV. This is essentially due to the measure of 
the contribution of the optical and infrared sources, which they 
find at the level of 28% of the total CXRB in this particular en- 
ergy interval (see Sect l4.1b . At odds with this measure, we do not 
find any significant feature in the X-ray spectrum of optically/IR 
detected sources, resulting in a estimate of the unresolved flux 
at 7 ke V which is four times higher than iLuo et al" I (|20TT1) 1 cr 
upper limit. 

In terms of the resolved fraction of the to tal CXRB, using th e 
Swift measurement of the total CXRB dMoretti et all |2009|) . 
our measurement of the unresolved components corresponds to 
-80% in the 2-10 keV band, with 95+3% at 2 keV and 70+16% 
at 7 keV (lower panel of Fig. [TOt . We note that this number 
should be kept with caution and does not really add any sig- 
nificant piece of news once the unresolved emission is directly 
probed and constrained. This is not onl y because the spectruni 
of the total CXRB i s highly uncertain ilRevmvtsev et all l2005t 
iFrontera et 311120071: iMoretti et ani2009h . but also because deep 
fields are limited to a small portion of the sky so they miss the 
rare and bright population that is responsible for ~20% of the 
total emission. 

To summarize, we can state that, compared to previous works 
and thanks to the low level of the Swift-XRT instrument back- 
ground in the hard band (>2 keV) we significantly reduced the 
uncertainties finding consistent results in terms of absolute flux 
measurement. After reducing the uncertainties our result left 
room for a small, but significant component still to be resolved. 



5.2. Comparison with AGN population modeis 

To check whether the unresolved spectrum can be accounted for 
by AGN population sy nthesis models, we compared our results 
with iGilli et all ( 120071) (G07) models 0, which provides the in- 
tegrated spectrum for different AGN subpopulations at given lu- 
minosity (L), redshift (z) and absorbing column density (Nh). 
To compare data and model, first of all, we need to assess which 
AGN subpopulation can be detected in the CDFS observation. 
Therefore, for a grid of six nn (lO^^-lO^^ cm"^), thirty -SIX L 
(lO'^'-lO'*^ erg s"' cm"^), and twelve z (0-6), we calculated the 
expected median flux, using the same spectral models as used in 
G07. For each bin flux we calculated t he detection proba bility by 
mea ns of the respons e function of the lXue et alj (1201 lb catalog. 
The iXue et all (1201 11) catalog is the result of the cross-correlation 
of three difi'erent catalogs built in three different energy bands, 
while the response function is only provided separately. We cal- 
culated the detection probability as the combination of the prob- 
ability in the soft and hard bands. To ease the calculation we ne- 
glected the possibility that a source can only be detected in the 
total band, while it is undetected either in the soft or in the hard. 
In our ROI this is only true for 12 out of 430 sources (2.7%). 
For example, considering AGN with photon index r= 1 .9 in the 
bin 10^2 <nH <10^\ 1.0<z< 1.5 and lO'*' '^ <L<IQ'^^-^ erg s"' 
cm"^, we expect fluxes of 6.0 x 10"'^ and 1.8 xlO"'^ erg s"' 
cm"^ in the 0.5-2.0, and 2.0-8.0 bands, respectively. The catalog 
response function (Table 7 and Fig. 23 of lXue et al.l(i201 II) ) gives 
the source detection probability at these fluxes in the two bands 
we are interested in, p.5o//=0.4 and p/,a,.j=0.2. We calculate that, 
for this AGN, the probability to be included in the CDF-S catalog 
is p^^,=1.0- (l-p.5o/,)x(l-p/,fl„/), that is 0.52. 
As in G07 a distribution of photon index was considered, we re- 
peated the same procedure using nine different slopes (between 
1.5 and 2.3) and weighting the final output by the same distribu- 
tion used by G07 that is a Gaussian centered on 1 .9 with standard 
deviation 0.2 . 

The G07 model is fu lly consistent (7% scatter) with the summed 
spectrum of all the in' Xue et"al1 (1201 ll) catalog (Fig.fTTTi. We also 
find that the same model accurately reproduces the soft part (< 
3 keV) of the spectrum of the unresolved emission, whereas it 
falls short in replicating the 3-7 keV emission, hinting that there 
are some missing hard sources. 

5.2.1 . AGN and galaxies 

To compare our data with the G07 model we ass umed that the 
emission of the 5683 optical/IR sources from the iGrazian et al.l 
(2006) catalog, subtracted from the unresolved spectrum, is 
dominated by non-nuclear sources, with a negligible contri- 



http : //www . bo . astro . it/~gilli/xrb . html 
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Fig. 11. In both panels the dark and light gray areas represent the spectrum of the unresolved and resolved CXRB, respectively. Left 
Panel: The yellow dashed and continuous lines show the expectations of the G07 model for the two components. Right Panel: 
The shaded area represents the total CXRB emission. Blue lines show the G07 model output with the CT (any z, any L, nH > 
24) contribution enhanced by a factor 2. Green lines show the G07 model with the only heavily CT contribution enhanced by a 
factor 5. Red lines show the case where the contribution of high z CT has been enhanced assuming a linear evolution with redshift. 
Continuous, long-dashed, and short-dashed lines are for CDFS unresolved, CDFS resolved, and total CXRB, respectively. 



bution from an AGN component. There is some evidence of 
this hypothesis. First, although we cannot measure the X-ray 
fluxes of the single sources (which are by definition fainter than 
the detection limit), we can estimate the mean luminosity of 
the sample, using the cataloged photometric redshifts. We find 
that a mean luminosity of 2x10-'^ erg s"' cm"^ in the 0.5-2.0 
keV band accounts for the measured total flux of the sample 
(7.25!f^^xl0"'^ erg s^' cm^^ deg^^ see Table [3]l. Such faint 
mean luminosity, together with the softness of the spectrum 
(photon index ~ 1.9), is a clear indication that the galactic emis- 
sion is prevalent with respect to the nuclear component. While 
the mean luminosity by itself can be an indication, it does not 
rule out the presence of some AGN. However, integrating (the 
extrapolation of) the X-ray galax y number-counts distribution 
measured bv iLehmer et al.l (l2012h in the 2.-10 keV band down 
to the flux (3x10"'' erg s~' cm"^) for which the source den- 
sity corresponds to the observed (~ 400,000 galaxies per deg^) 
and accounting for the catalog completeness function, we find 
a value of 6.5 10"'^ erg s"' cm"^ deg"^, which is very close to 
and consistent with our measure. Moreover, a completely inde- 
pendent estim ate of the galaxy X-ray flux density comes from 
iDijkstra et al. ^('2012) who estimate the contribution of the galax- 
ies to the CXRB. Lacking an accurate measure of the mean 
galaxy spectral slope in the hard band, they calculated this con- 
tribution by assuming different mean photon indices and using 
the star formation history of Hopkins & Beacom (2006) together 
with the star formation rate to X-ray luminosity conversion of 
iMineo et al.l (j2012). This estimate is particularly suitable for our 
purposes as we have a direct measure of the galaxy mean spec- 
tral slope. In the case of r=1.9, they estimated that the galaxy 
density flux is 1.2+0.3x10"'^ erg s"' cm^^deg"- in the 2-10 
keV band, which is slightly higher than our measure and per- 
fectly consistent with our assumption. These arguments suggest 
that the emission of the X-ray undetected optical/IR sources of 
iGrazian et al. I (Hop) catalog is mostly of galactic origin with a 
negligible nuclear contribution. 



5.3. Some speculations 

We examine the possibility that the discrepancy between the un- 
resolved spectrum and the G07 model could be explained by 
the cumulative emission of a number of Compton-thick (CT) 
which are not included in the model. From an observational 
point of view the CT lu minosity function is highly uncertain 



ay 

f 



both in the local Universe (Delia Ceca et alJl200 8: Treist er et al 
'2009';'Banantvne et al."201 iVSevergnini et al.'2012; Aiello et al 
2012) and at higher redshift (Tozzi et al. 2006; Daddi et al. 20G 
Alexander et al. 2008; Fiore et alJl2009t lAlexander et al.ll20r 
In G07 the CT luminosity functions at different redshift are in- 
directly calculated from the difference between the total CXRB 
and the integrated emission of Compton-thin AGN (21<Nh < 
24), under the assumption that the number of mildly CT objects 
(log Nh - 24.5) is equal to that of heavily CT objects (log Nh 
- 25.5) and that they have the same cosmological evolution of 
Compton-thin AGN. 

While the total integrated emission of the CT population by def- 
inition agrees with the total CXRB, as we said in the previ- 
ous section, we found some discrepancies with the constraints 
placed by the unresolved emission spectrum. Multiplying the 
whole CT contribution (at any redshift and luminosity) by a fac- 
tor 2.0 would reconcile model and the observation in the 1.5-7 
keV band, but, at the same time, would break through the spec- 
trum of the total CXRB at harder energies ~ 20-30 keV, where 
the CT contribution is less affected by absorption (blue line in 
the right panel of Fig. [TTT i. 

Enhancing the contribution of the only heavily CT population 
by a factor five significantly reduces the distance between data 
and model; but the number of these sources, which is required 
to make the model consistent with the data under 7 keV, would 
result in over predicting the total CXRB at higher energy (green 
line in right panel of Fig. fTTI) . 

On the other hand, high-redshift (z >2) CT are expected to peak 
at lower energies and to supply less flux at energies > 20 keV. 
Indeed, while in G07 model a constant fraction of obscured AGN 
with redshift is assumed, we find that, hypothesizing a linear 
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Fig. 12. The comoving space density of the CT AGN. We report 
different results for different redshift and luminosity intervals. 
Different colors indicate different luminosity bins. Continuous 
lines show the output of the GOV model; dashed lines represent 
the same outputs modified by applying a linear positive evolu- 
tion in redshift. The green dotted line displays the output of the 
composite model by , Draper & Ballantvne (2010) 



positive evolution with redshift (1+z) in the total number of CT 
objects, corresponding to a (l+z)" with a~ 0.5 in terms of CT 
fraction, would solve the discrepancy between model and ob- 
servation. With this assumption the overall AGN contribution is 
still consistent with the CXRB constrain at higher energies (red 
line in the right panel of Fig. [TTT i. 

In Fig. [12] we report different measurements of the CT space 
density performed with different techniques (most in the infrared 
band) at different redshift, for different X-ray luminosity inter- 
vals, together with the output of the G07 model and imposing 
a linear evolution as suggested by our observation. This shows 
that such an evolution of the absolute number of CT AGN is still 
consistent with the large scatter of the observation results. 
Our solution would require an increase in the CT fraction with 
redshift. The change in the obscured AGN fraction with red- 
shift is a debated issue. A positive evolution has been ob- 
served in some studies jL a Franca et al. 20051 iTreister et al.l 
120061: iBallantvne et al.l 20Q6; Hasineer 2008), whereas ot hers do 
not fi nd any significant variation (Ueda et al. 2003; Akvl as et al.l 
l2006l) . While these studies adopted a value of ~10^'cm"^ as 
pivotal column density to separate obscured and unobscured 
sources, our data seem to indicate a stronger evolution for only 

the CT population (>102''cm"2)^ 

A similar result was found by iBrightman & Uedal (1201 2h that 
shows a significant increase in the only CT fraction, from 20% 
in the local universe up to 40% at z-1-4. An higher number 
of hig h-redshift CT have recently been claimed by Gilli et al. 
(120111) who find a CT AGN at z=4.76 in the CDF-S area, the 
detection probability being less than one. A strong positive evo- 
lution of the only C T population, has been also discussed in 
ITreister et al] (1200 9*) as a possible explanation to reconcile the 
difference between IR observations and the expectations from 
AGN population models normalized to reproduce the local CT 
density. 



From the viewpoint of the population m odeling, evolution of 
the only CT population has been found bv lDraper & B allantvnd 
(2010) as a natural consequence of describing the CT AGN by 
means of a physically motivated Eddington ratio distribution: CT 
AGN are a composite population accreting either at >90% or 
<1%. Indeed in this model the predicted evolution of the num- 
ber density of CT AGN seems to follow the one we propose, 
although the normalization is a factor 2 lower (Fig.fTSI). 
In the framework of AGN-galaxy coevolution, a positive red- 
shift evolution of heavily obscured AGN is in line with the idea 
that in the high-redshift Universe SMBH are hosted by gas-rich 
galaxies. In particular, one of the scenario proposed for the for- 
mation and evolution of SMBH and galaxy assumes that major 
merger can trigger star formation and black hole growth at the 
same time. Depending on the kind of interaction, different ob- 
scuration geometries could be possible around the black hole, 
but in any case this would increase the probability of observing 
a SMBH through a very large amount of gas and dust. 

6. Conclusions 

We exploited the low and predictable instrument background of 
the Swift XRT telescope to perform the spectroscopy of the un- 
resolved X-ray emission in the CDF-S. We found a faint, but 
significant unresolved component that can be modeled by a very 
hard power-law with photon index r=0.1±0.7 and a flux den- 
sity of 5 xlO"'^ erg s"' cm"^ deg"^ in the 2.0-10. keV band, 
coiTesponding to 20% of the total CXRB. With respect to pre- 
vious works we significantly improved the accuracy over the 
1.5-7.0 keV band. Our measure is in very good agreement with 
what is expected by the G07 AGN population model in the 1.5- 
3.0 keV. In the hard band (3-7 keV) the same model falls short 
when replicating the observed spectrum, pointing toward some 
missing very hard sources. This discrepancy can be solved hy- 
pothesizing a positive evolution with redshift of the contribution 
Compton-thick AGN population. 
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